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ABSTRACT 



Context. Radio, X-ray, and gamma-ray observations provide us with strong evidence of particle acceleration to multi-TeV energies 
in various astrophysical sources. Diffusive shock acceleration is one of the most successful models explaining the presence of such 
high-energy particles. 

Aims. We discuss the impact of inverse Compton losses on the shock acceleration of electrons that takes place in radiation dominated 
environments, i.e. in regions where the radiation energy density exceeds that of the magnetic field. 

Methods. We perform a numerical calculation, including an energy-loss term in the transport equation of accelerated particles. 
Results. We discuss the implications of this effect on the hard X-ray synchrotron and gamma-ray inverse Compton radiation, produced 
by shock-accelerated electrons in young supernova remnants in the presence of large radiation fields (e.g. in the Galactic centre). We 
also discuss possible implications of our results for clusters of galaxies and gamma-ray binaries. 

Conclusions. We demonstrate that the inverse Compton losses of electrons, in the Klein-Nishina regime, lead to spectra of ultra- 
relativistic electrons that may significantly differ from classical diffusive shock acceleration solution. The most prominent feature is 
the appearance of a pile-up in the spectrum around the cut-off energy. 

Key words, shock acceleration - radiation mechanisms: non-thermal 
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1. Introduction 

Diffusive shock acceleration (DSA) is the most successful and 
widely accepted mechanism for explaining the acceleration of 
cosmic rays in many astrophysical environments. Despite the ba- 
sic physic s of DSA being very robust and reasonably well under- 
stood (see lBlandford & EichleJl987l and lMalkov & Drurvl200ll 
for a review), some crucial issues still need to be addressed. 

In its simplest, test-particle approach the theory of DSA pre- 
dicts an energy spectrum for the accelerated particles, which is 
a featureless power law f{p) oc p^" up to a maximum energy 
where the spectrum cuts off. For strong, nonrelativistic shocks, 
where the mechanism is believed to operate very efficiently, the 
slope of the spectrum below the cut-off converges to the canoni- 
cal slope a - A, independently of the details of the acceleration. 

The value of the maximum energy attained by particles is de- 
termined by the competition between the acceleration time and 
the shortest of the three time scales: /) particle escape time from 
the accelerator, //) particle energy loss time, and ///) age of the 
accelerator. Therefore, while the energy spectrum below the cut- 
off is almost universal, the spectral shape of the cut-off depends 
dramatically on both the details of the acceleration mechanism 
and the physical process determining the particle maximum en- 
ergy. In the case of electrons, which we consider in this paper, 
the maximum energy is in most cases limited by radiative syn- 
chrotron and inverse Compton losses in the ambient magnetic 
and photon fields. 

The missing piece of information in the model is the dif- 
fusion coefficient of particles close to the shock, which deter- 
mines the acceleration rate. For a given diffusion coefficient, 
the problem is well-defined and can be solved once the value 
of the magnetic field and the spectrum of the ambient radiation 
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are specified. Although there is some uncertainty about this co- 
efficient, one can reasonably assume that, due to the high level 
of turbulence expected close to the shock, the mean free path 
of particles is close to their Larmor radius r^. Under these cir- 
cumstances, diffusion proceeds close to the slowest possible (so- 
called Bohm) rate with diffusion coefficient D ~ rz,c/3. 

Several approaches can be found in the literature to the prob- 
lem of elec tron acceleration at shocks in the presence of radia- 
tive losses. Bulanov & Dogiel ('19791). IWebbet all ( 11984 . and 
Heavens & Meisenheimer (1987) considered the case of a dif- 
fusion coefficient constan t in m omentum and space. Recently 
IZirakashvili & AharonianI (12007!) numerically solved the prob- 
lem for D(p) with an arbitrary momentum dependence. All these 
approaches are limited to the case of dominant synchrotron 
losses. 

In this paper, we numerically solve the problem for a general 
form of both the diffusion coefficient and the energy-loss rate. 
In particular, we focus on the case, never considered before, in 
which the accelerator is embedded in a strong radiation field, 
characterised by a much higher energy density than the magnetic 
one. Under this circumstance, inverse Compton losses dominate 
synchrotron losses and, at the highest energies, modifications 
to the electron spectrum are expected with respect to the case 
of synchrotron/Thomson losses, due to the transition between 
Thomson and Klein-Nishina regimes. This, in turn, strongly af- 
fects the spectrum of the radiation emitted by accelerated parti- 
cles. We calculate the exact shape of the spectrum in the whole 
energy range and show that, in the case of Klein-Nishina losses, 
the particle distribution at the shock has a broad cut-off, due to 
the shallow energy dependence of the loss rate. Another impor- 
tant effect of Klein-Nishina losses is to harden the downstream 
electron spectrum close to the maximum energy, leading to the 
formation of a pronounced pile-up. Such a feature can be ob- 
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served in the photon spectrum, in particular for the synchrotron 
emission. The effect on the inverse Compton emission is less 
pronounced because the Klein-Nishina cross-section intervenes 
twice, in opposite directions, to harden the electron spectrum and 
to soften the photon one leading to an almost exact compensa- 
tion. 

We describe the details of our model in Sect. |2] and dis- 
cuss the main characteristics of the calculation by applying it 
to an ideal example in Sect. |3] In Sect. |4] we then apply our re- 
sults to several astrophysical environments where the radiation 
field-energy density dominates, namely to supernova remnants 
located close to the Galactic centre, shocks in binary systems 
and accretion shocks around massive clusters of galaxies. For 
these systems we evaluate the spectrum of the accelerated par- 
ticles, together with the spectra of both the emitted synchrotron 
and inverse Compton radiation. We summarise in Sect. |5] The 
presence of a strong pile-up at the high-energy end of the syn- 
chrotron spectrum is the most remarkable feature. 

2. The model 

We consider a plane-parallel shock where the fluid moves (in the 
shock rest-frame) along the x-axis from -oo far upstream, to -i-oo 
far downstream, and the shock is located at x = 0. We assume 
the velocity of the shock to be nonrelativistic. 

The transport equation for the particles distribution function, 
including the presence of energy losses, reads as 



df(x, P^t) df(x, p,t) d 1^^ ^^^ df{x,p,t) \ pdu df(x,p,t) 
dt dx dx\ ' dx j 3 dx dp 



2.1. Electron spectra 

We solve Eq. ([1) using a finite difference numerical scheme, im- 
plicit in time to guarantee unconditional stability and explicit in 
momentum. To solve the problem we need to insert the bound- 
ary conditions at the shock and at upstream/downstream infinity. 
For the boundary at the shock we consider again Eq. ([T]) and in- 
tegrate it between 0- immediately upstream and 0+ immediately 
downstream, obtaining 



-(Ml - U2)p— = D2 — \ - Di—\ + Qo6(p-po) 
3 op ox 12 ox II 



p^ dp 



— —(P'-(^^PWx,p,t)) 



Q(x,p), 



(1) 



where L(x,p) - -p is the loss rate taken to be positive and 
Q{x,p) is the injection term; u represents the bulk velocity of 
the plasma in the shock rest-frame. In the following we as- 
sume that the diffiision coefficient in the up and downstream re- 
gions Di 2 does not depend on x and has the functional form 
Dip) - Dqi/ (although any expression for Di 2{x, p) can be eas- 
ily implemented in the code) and that injection happens at the 
shock surface as a delta function in energy: 

Q(x,p) = Qa6(x)6(p - pq). 

The formulation is general so that it can be applied to both 
protons or electrons, plugging in the relevant energy-loss chan- 
nels, namely proton-proton and proton-gamma interactions, for 
the hadronic channel, and synchrotron and inverse Compton (IC) 
emission for the leptonic one. In this paper we study the acceler- 
ation of electrons undergoing synchrotron and inverse Compton 
losses, in the case where the IC is dominant. This is the case 
if acceleration takes place in a region where the energy den- 
sity of the ambient radiation field exceeds the magnetic one: 
Urad » l{%n). In particular we are interested in studying the 
effect produced on both the electron and the photon distributions 
when IC losses proceed in the Klein-Nishina regime, namely 
when eEI(mc^)^ > 1 with E being the electron energy, m the 
electron mass, and e the ambient photon's energy. 

In the following we assume that the background radiation 
field is isotropic and that the background magnetic field is con- 
stant in the upstream and downstream regions and the two are 
related by B2 - ^B\, where ^ is the compression factor Such a 
parameter varies between two values: ^ = 1 for a pure parallel 
shock and ^ - A- for a perpendicular shock. 



(2) 



where the subscript 1 refers to the upstream region and 2 to the 
downstream one. At +00 we set fix, p) - 0. 

It is convenient to introduce two dimensionless variables in 
place of p and x. The first variable is p/p*, where p* is a parame- 
ter that estimates the cut-off momentum for the electron distribu- 
tion, evaluated by imposing equilibrium between the momentum 
gain in one acceleration cycle. 



^Pacc = 



4p (mi - U2) 



3 C 

and the momentum loss per cycle, 
Apioss = Liip)Ati + L2ip)At2, 



(3) 



(4) 



where the mean residence ti me in the up and do wnstream regions 
is Afi,2 = 4Di,2(p)/cMi,2 (see lWebbetal]|l984l) . 

For the spatial coordinate we operate a change of variable 
both up and downstream. To reduce the wide range in |x| from 
to 00 to one more numerically feasible, we use an exponential 
variable: 



exp 



_|x| 
xo 



(5) 



where xq is a characteristic length scale. The new variable now 
conveniently ranges from to 1 . In the upstream region, except 
for energies close to the cut-off, the accelerated particles can 
propagate one diffusion length ahead of the shock. The situa- 
tion is different downstream, where the maximum distance from 
the shock a particle can reach, in steady state, is determined by 
energy losses. Thus, a reasonable choice for the length scale xo 
is 



and 



Diip) 

Ml 



X0,2 = U2TL, 



(6) 



(7) 



for upstream and downstream, respectively. The energy loss time 
is given by tl - pi Lip). 

Once the solution /(x, p) of Eq. ^ is found, we evaluate 
the total spectrum integrated over space, Fip) - J fix, p)dx, 
which is the quantity required to evaluate the radiation spec- 
trum emitted in the shock region. It is worth mentioning that the 
steady state integrated spectrum, below the cut-off energy, fol- 
lows a power law. In the downstream region, in case of Thomson 
cooling, the differential spectrum in momentum space has index 
6 - a + I =5. Upstream, the maximum distance from the shock 
a particle can travel is essentially its diffusion length; therefore, 
for Bohm-like diffusion, this results in 6 = a - 1 = 3. 
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2.2. Photon spectra 

Once the electron spectrum is obtained, we can calculate the 
spectra of the radiation emitted both via synchrotron and IC. For 
synchrotron radiation the energy flux at an energy e is given by 



<D(e) 



hmc^ 



(8) 



where Cc - hvc is the energy corresponding to the critical fre- 
quency Vc - lieBp^ l(Annr'c^) and Kiejec) is the emission pro- 
duced by the single electron of momentum p, charge e, and mass 
m; h represents the Planck constant. 

In our case F{p) is the total electron spectrum integrated up 
and downstream. The exact expression for the kernel function 
K (e/ec) in the case of a turbulent magnetic field was derived 
in lCrusius & Schlickeiseil ( [T9 86). With several percent accuracy, 
th is can be approximated by the an alytical expression obtained 
in lZirakashvili & Aharonianl ( l2007h : 



V(e/e,)-2/3 + (3.62/71)2 



(9) 



We adopt this simpler expression. 

The inverse Compton energy flux for an isotropic distribu- 
tion of soft photons «(e') upscattere d by a pop ulation of elec- 
trons with spectrum F{p) is ( Blume nthal & Goul d 1970.) : 

, Ine'^e C , , C n(e')de' 
0(e) = J dpFip) J 



1 



2qlnq + {l + lq){\ - q) + - 



2 pc(pc - e) 



(1-?) 



(10) 



where 



q = 



p-^{pc - e) 



For simplicity, in the following we assume that the back- 
ground radiation field is a black body, though these calculations 
may easily be extended to any background radiation spectrum. 
Depending on the environment we want to model, we consider 
either a single or a superposition of multiple Planck distribu- 
tions, introducing a dilution factor r\ for each; 



«(e') = 7? 



1 



(11) 



where k is the Boltzmann constant and T the black body temper- 
ature. 



3. Results 

We present here the results of our calculation by means of an 
ideal case to illustrate the features and characteristics of the 
problem we are studying. We consider a nonrelativistic shock 
expanding in a medium where the magnetic field value is = 
10 plG. For such a field, the magnetic energy density is l'S,n - 
2.5 eV/cm-^. To keep a limited number of model parameters, in 
this example we assume the same value for the magnetic field 
upstream and downstream (i.e. ^ = 1 , strictly valid for a perfectly 
parallel shock). In the next section, when we apply our results to 
astrophysical objects, we include the compression of the mag- 
netic field at the shock. We assume that the accelerator is em- 
bedded in an isotropic background radiation field. We further as- 
sume a diluted black body radiation spectrum (i.e. with a reduced 
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Fig. 1. Energy-loss rates and shock acceleration rates for a mag- 
netic field B = 10 yuG and a diluted black body radiation field at 
a temperature T - 3500 K. All curves have been multiplied by 
£■ 2 so that Thomson losses correspond to horizontal lines. Thick 
lines: synchrotron loss rate (dot-dashed), inverse Compton loss 
rate (dashed) and the sum of the two (solid). Thin curves: ac- 
celeration rates for ui - 1000 km/s (dashed), u\ - 3000 km/s 
(solid) and ui - 10000 km/s (dotted). 



energy density) at a temperature of kT - 0.3 eV (i.e. in the near 
infrared energy band), corresponding Xo T - 3500 K. We then 
assume the radiation energy density to be Urad - 5x10^ eV/cm-', 
a factor 2x10"^ greater than the magnetic one. (As we point out 
in the next section, this high value is simil ar to the one mea sured 
in the < 1 pc Galactic centre region (Davi dson et al.|[T992l) ). 

In Fig. [T] we plot the energy-loss rates due to synchrotron 
losses in a magnetic field of 10 pG and inverse Compton scatter- 
ing off" a radiation field at 3500 K, as a function of the particle 
energy. The sum of the two is also represented, as well as the 
acceleration rates for DSA in the same magnetic field and for 
three different values of the shock velocity {u\ - 1000 km/s, 
M] - 3000 km/s, and u\ - 10000 km/s). In this figure all curves 
have been multiplied by E^^ so that Thomson losses correspond 
to horizontal lines. The intersection between the acceleration 
rate curve and the total energy loss rate indicates the point in 
the particle energy spectrum where the acceleration due to DSA 
is compensated by radiative losses so the cut-off sets in. 

Figure[T]shows how the energy dependence of energy losses 
changes behaviour At low energies (E <s: {mc^)^/epi,), inverse 
Compton proceeds in the Thomson regime (dE/dt oc E^), but at 
high energies the losses enter the Klein-Nishina regime and the 
dependence of the process on energy changes (dE/dt ~ InE). At 
even higher energies synchrotron losses become dominant and 
an E^ law for the loss rate is recovered. 

We solve the transport equation for the accelerated electrons 
choosing the value of 3000 km/ s for the shock velocity. For such 
a choice of parameters, losses proceed deep in the KN regime. 
The shock is assumed to be strong, i.e. compression ratio = 4, 
and we also assume a Bohm-type diffusion coeflicient, namely 
D - D^)p^ with P - I and Do = c^/{3eB). In this example we 
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Fig. 2. a) Electron spectrum at the shock location in the case of 
inverse Compton dominated losses (solid line) and synchrotron 
cooling dominated case (dashed line) for a magnetic field B - 
10 /iG equal up and downstream, shock velocity mi - 3000 km/s 
and compression ratio R = 4. The background radiation field is 
assumed to be a diluted black body at T = 3500 K with en- 
ergy density Urad = 5 x 10"* eV/cnr'. The momentum scale for 
each curve is normalised to the parameter p*. The normalisation 
on the y-axis is in arbitrary units, b) Electron distributions inte- 
grated over space. In the case of the IC dominated case (solid 
line), we also show separately the two components upstream 
(dotted) and downstream (dash-dotted). The dashed line repre- 
sents the synchrotron dominated case. 



make the assumption that the radiative energy loss time scale is 
much shorter than any other loss time scale in the system that 
could limit the acceleration, such as escape or the finite lifetime 
of the accelerator, at all the energies we consider. Therefore we 
look for the steady state spectrum. The resulting electron energy 
spectrum at the shock location is shown in Fig.|2^, together with 
the spectrum expected for the same values of the parameters but 
suppressing the radiation energy density (i.e. when synchrotron 
losses dominate). For an easier comparison we rescaled the mo- 
mentum scale for each curve to its cut-off value p*. Of course the 
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Fig. 3. a) Inverse Compton photon distribution produced by the 
electron spectra in Fig.|2]3 on a diluted black body atT - 3500 K 
for the case of dominant IC losses (solid line) and for the case 
of dominant synchrotron cooling (dashed line), b) Synchrotron 
emission from the two electron distributions in Fig.l^J, the mag- 
netic field value is 10 fiG. 



actual value of p* is different in the two cases. The normalisa- 
tion on the y-axis is in arbitrary units. The shape of the spectrum 
is modified by the decrease of the IC cross-section at high en- 
ergies: the cut-off region is broader and the decay is shallower, 
compared to the synchrotron case. The result may be understood 
most easily from the rates in Fig. [1] which highlights the weak 
dependence on energy of Klein-Nishina losses. 

To obtain the overall radiation produced by this population 
of electrons we have to consider the spatially integrated spectra 

Fi - f(x, p)dx and F2 - f{x,p)dx. Figure |2]3 shows 
the total spectrum F{p) - F\{p) + Flip) resulting from the inte- 
gration for both cases of IC dominant losses and of synchrotron 
dominant losses. As can be seen, the most pronounced feature 
appears in the integrated spectrum: the softer dependence on 
energy of the energy loss rate in the KN regime, as compared 
to the Thomson regime, causes a significant pile-up around the 



G. Vannoni et al.: Diifusive shock acceleration in radiation dominated environments 



5 



cut-ofF energy. For this case, we also plotted the upstream and 
downstream distributions separately to show that the major con- 
tribution to the total spectrum is provided by the downstream 
distribution. In this region all particles are eventually advected 
away from the shock with the same velocity U2- They also dif- 
fuse to a distance that depends on their energy and on time as 
xd - '\J4-D(p)t. For particles with energy below the cut-ofF en- 
ergy, advection is always dominant and the maximum distance 
to which they can propagate is determined by their loss time 
scale and is Xad - M21"l- The diffusion length at the same time is 
given by xo = yjAD{p)Ti. Both the length scales increase with 
tl- In particular, the two are comparable at the cut-ofF energy 
where = ^ Dju^. At lower energies diffusion is negli- 
gible and advection sets the maximum propagation distance. In 
the first approximation, the particle integrated spectrum down- 
stream can be written as F2{p) - foxi, with xi - mm{xaii,xi,]- 
For Thomson scattering ti oc £""', determining that the higher 
the particle energy, the smaller the distance it can travel away 
from the shock before losing its energy. For this reason the in- 
tegrated spectrum in momentum space downstream, at energies 
below the cut-oiF where xi - Xad, is F2ip) = /quiTl p^"^^ ■ In 
the case of KN losses we have an opposite trend: ti oc E/lnE, so 
that high-energy particles can propagate farther than low-energy 
ones, and the resulting spectrum is harder This effect produces 
the pile-up when integrating the spectrum over space. In the case 
we are considering in this section, we assume that the loss time 
scale is the limiting one at every energy in the range we consider 
In real physical environments this may not be the case. At low 
energies, where the Thomson regime is recovered also in the IC 
case, the time scale of energy losses increases with decreasing 
energy, therefore it could at some point become more than the 
age of the system (Tags) itself. The maximum energy of the ac- 
celerated particles is still determined by the energy losses, but 
the low-energy part of the spectrum does not have enough time 
to cool and the spectrum results in F2 - foU2Tage ^ P 

With the total spatially integrated spectrum obtained, we cal- 
culate the radiation emitted. Our results are shown in Fig. [3] 
Panel a shows the inverse Compton radiation resulting from 
the upscattering of the background radiation, when losses are 
Compton-dominated, compared to the dominant synchrotron 
losses case. The spectrum pile-up is not so remarkable in the 
first case because the efFect of the Klein-Nishina cross-section 
acts twice: while hardening the electron's distribution, it softens 
the IC photons' one, so that the two effects almost cancel each 
other out. 

On the other hand, the accelerated particles' features have 
a sharp imprint on the synchrotron spectrum (Fig. which 
carries direct information on the electron distribution. Compared 
to the case of pure synchrotron cooling, the emission due to an 
electron distribution shaped by IC losses presents a pile-up of 
three orders of magnitude around the cut-ofF, which appears at 
keV energies for the chosen parameters. This characteristic is a 
remarkable signature to distinguish the two different scenarios 
for electron radiative losses. 
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Fig. 4. a) Energy loss rates and acceleration rates, multiplied 
by E^^, as a function of the electron energy E, for a super- 
nova remnant in the Galaxy. The synchrotron rate (thick dot- 
dashed line) corresponds to a magnetic field Z? = 10 fiG, the 
IC curves are calculated for three black bodies at temperature 
T = 2900 K (optical), thick long-dashed line; T = 29 K (far in- 
frared), thick short-dashed line; T = 2.7 K (CMB), thick dotted 
line. The thick solid curve represents the sum of the four loss 
contributions. Thin lines: acceleration rates for Z? = 10 /iG and 
shock velocity of 1000 km/s (dotted), 3000 km/s (solid), and 
10000 km/s (dashed), b) Same as in panel a but for a magnetic 
field Bi - 100 fiG upstream and B2 - 4-Bi downstream, and 
three black bodies at T = 35000 K (UV/Optical), T = 3500 K 
(NIR), and T = 70 K (FIR). 



4. Applications 

The results we obtained in the previous section may have broad 
astrophysical applications. To demonstrate this, we apply the 
results of our calculations to three specific astrophysical envi- 
ronments where strong nonrelativistic shocks may form and in 
which the radiation energy density may dominate over the mag- 
netic field energy density. The three cases we consider are: /) a 



supernova remnant in the Galactic centre region, //) the accretion 
shock surrounding clusters of galaxies and ///) a shock in a mi- 
croquasar jet. For each of the considered cases, we evaluate both 
the electron spectrum and the spectrum of the emitted radiation. 
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Fig. 5. a) Electron spectrum at the shock for a velocity of 
3000 km/s and a magnetic field Z?i = 100 /iG. The solid line 
is obtained in the case of dominant IC losses with the same pa- 
rameters as in Fig. the dashed line corresponds to the case 
when the radiation energy density is suppressed and synchrotron 
cooling dominates, b) Integrated spectra for the same values as 
the parameters of panel a. The three curves correspond to differ- 
ent maximum times of acceleration, as reported. 



4.1. SNR in the Galactic centre 

Several observations in the radio and X-ray bands have con- 
firmed supernova remnants (SNR) shock s as powerful accel - 
erators of electrons up to tens of TeV (iBamba et al.l l2003h . 
We first consider a SNR in the galactic disc. The radiation 
field in the Galaxy consists of three distinct components: the 
optical/near-infrared (NIR), the far infrared (FIR) and the cosmic 
microwave background (CMB) radiation, with an energy den- 
sity o f ~ 0.5, ~ 0.1, ~ 0.25 eV/cm^, respectively (iMathis et al.1 
Il983h . In our calculation we evaluated an effective tempera- 
ture for these fields from their peak energy and then assumed 
they can be approximated with a diluted black body distribution. 
The Galactic magnetic field has a measured value of ~ 3 fiG 
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radiation spectra. 



(see IWidrowll2002l for a review), which corresponds to an en- 
ergy density of ~ 0.23 eV/cm^. This implies that the aver- 
age magnetic-field energy density is comparable to the radiation 
one. Moreover, the magnetic field can be significantly amplified 
in the presence of a shock that is efficiently accelerating par- 
ticles (iBell & Lucekll2001h . so that, in the acceleration region, 
t^mog >> Vrad results. Therefore synchrotron losses dominate 
the inverse Compton (Fig. |4^). In this case the contribution of 
the IC to the energy-loss rate is negligible and the cut-off energy 
and shape are determined by synchrotron cooling. 

Things are dramatically changed, though, if we consider an 
SNR in the Galactic centre region, where both the magnetic and 
the radiation fields are much stronger than the average Galactic 
values. Analogous to what is done in iHinton & AharonianI 
(2007), we consider the inner 1 pc of the Galaxy, where the 
values of the radiation energy density are 5 x 10"* eV/cm^ at 
kT = 3.0 eV (UV) and kT = 0.3 eV (NIR), and 5 x 10^ eV/cm^ 
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at /tr = 6 X 10"^ eV (FIR) (iDavidson et al.lll992h . To have an 
equal energy density in the magnetic field, it is required that 
B ^ 500 fiG, any lower value implies the dominance of IC 
(KN regime) losses in determining the cut-off in the electron 
spectrum. We choose Z?i = 100 fiG, a value of the magnetic 
field compression ^ = 4, appropriate to Alfvenic turbulence, 
and a Bohm-type diffusion coefficient. For a completely ran- 
domised magnetic field, the compression factor would approach 
^ ^ VTT. As shown in the Fig.|4}3, for a velocity of 3000 km/s, 
Bi = 100 yuG and B2 = 4Bi, the cut-off energy is 10 TeV. 

As discussed in Sect. |3] depending on the age of the system, 
the particles in the low-energy part of the interval we consider 
may not have enough time to cool. Keeping the age of the histor- 
ical SNR as a reference, which are known to be efficient electron 
accelerators, we performed our time-dependent calculation and 
present here the results for T^ge = 100 yr, 500 yr and finally, as a 
comparison, the steady state case. 

At the shock location, steady state is reached quickly so that 
the three cases overlap perfectly. The result is shown in Fig. |5^ 
where we compare it with the one obtained for pure synchrotron 
losses, with the same values of the other parameters. The fea- 
tures are those pointed out in Sect. [3] the cut-off region is broader 
and the decay is not as steep as for synchrotron cooled electrons. 

In Fig.|5j) we plot the spatially integrated spectra for the three 
different time scales chosen. The up and downstream contribu- 
tions are also plotted for the steady state solution to show the na- 
ture of the substructure present around the cut-off energy. In the 
integrated distribution, the modification of the spectrum due to 
KN effects can be dramatic. In our example the hardening of the 
spectrum right below the cut-off causes a pile-up of more than 
one order of magnitude, compared to the solution for Thomson 
losses, which would correspond to a horizontal line in the plot. 
The shape of the spectrum around the cut-off is determined by 
losses, therefore the same as for the three curves. At lower en- 
ergies, the point where losses set in depends on the time. For 
those energies where ti > T^ge, we obtain the uncooled spec- 
trum F(p) cc p^*. We plot the steady state solution as a test to 
show that our method recovers the right solution oc p^^ at low 
energies where losses proceed in the Thomson regime. 

Once we have the electron spectrum, we can calculate the 
photon distribution. Our results for IC are shown in Fig. |6^, 
where we plot the spectra for the different time scales we are 
considering. Each of the curves represents the sum of the con- 
tributions of the three components of the seed field (FIR, NIR, 
UV). For the steady state solution, the overall radiation spec- 
trum shows an almost flat profile below the cut-off due to the 
compensation of the "double action" of the Klein-Nishina ef- 
fect. The other two curves harden in the low-energy part due to 
the slope of the electron spectrum at low energies. In addition, 
we see from Fig. |4j3 that the contribution of the three radiation 
fields dominate in different energy ranges. The higher the tem- 
perature of the black body, the lower the energy of the electrons 
required to enter the KN regime. For this reason, at early stages 
of the evolution (100 yr, dotted curve in Fig.|6h) when the low 
energy part of the electron spectrum is not yet contributing to 
the photon emission, the contribution of the two high tempera- 
ture black bodies at 3500 K and 35000 K is almost zero and the 
total IC radiation is lower than at later times. 

Fig lSti shows the synchrotron emission. Here the situation is 
different: as expected the radiation spectrum reproduces the fea- 
tures of the electron one. A pronounced pile-up appears around 
the cut-off energy for all the three cases. 
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Fig. 7. Thick lines: energy loss rates (synchrotron plus IC) for 
typical conditions in a galaxy cluster accretion shock. The back- 
ground radiation field is provided by the CMB and the val- 
ues of the magnetic field upstream are reported (downstream 
B2 - 4Z?i). Thin lines: acceleration rates for a shock velocity 
of 2000 km/ s and the line type refers to the corresponding value 
of Bi. 



In this treatment we have assumed a value Bi = 100 juG, 
reasonable for the Galactic centre region, although we should 
note that the magnetic field in this region is quite uncertain. The 
effect of a possible amplific ation of the magnetic field, accord- 
ing to iBell & LuceS jioOll) and .Bell C2004) . may lead to val- 
ues of the magnet ic field upstream of sev eral hundreds of micro- 
Gauss and more (* Uchivama et al.ll2007l) . Moreover, in the case 
of a modification of the shock due to the presence of accelerated 
particles, the total compression ratio can exceed the value of 4, 
therefore leading to an additional increase in the magnetic field 
value downstream. As a check, we performed our calculation 
for a compression ratio R -1 and a magnetic field compression 
^ - 1 and found that the effect of IC losses ceases to be sig- 
nificant, as compared to synchrotron cooling, for values of the 
upstream magnetic field exceeding 300 jiG. 

4.2. Clusters of galaxies 

Another example of a radiation-dominated environment where 
a shock can form is represented by clusters of galaxies. Rich 
clusters of galaxies are the biggest virialised structures in the 
Universe, with typica l sizes of a fe w Mpc and masses up to 
IO'^Mq or more (see I Sarazinlll988l for a review). An expand- 
ing shock wave, called the accretion shock, forms at the clus- 
ter boundary and carries the information of virialisation out- 
ward (Bertsch inger 19 85). The infalling matter crosses the shock 
surface at a speed roughly comparable to the free fall velocity, 
namely. 
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dominant synchrotron cooling, b) Integrated spectra for the IC 
dominated case. The upstream and downstream contributions are 
shown. 



so that the shock velocities are comparable to the one found 
in supernova remnants. Parti cle acceleration i s expected to hap- 
pen at accretion shocks (see Blasi et al."2007' for a review) with 
consequent emission of inverse Compton gamma rays, possi- 
bly detectable in the near future by GLAS T or ground-based 
Cherenkov telescopes (iGabici & Blasi"2004l). 

The detection of a tenuous and diffuse synchrotron radio 
emission from about one third of the rich clusters of galaxies 
dFeretti & JohnstoiT Hollittll2004HCarilli & Tavlo3l2002h reveals 
the presence of a magnetic field of a few /zG in the intraclus- 
ter medium. The value of the magnetic field in the cluster out- 
skirts, where the accretion shock propagates, is unknown, due 
to lack of radio measurements, but it is likely to be smaller or 
at most comparable to the one measured in the inner regions. 
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Thus, the magnetic-field energy density at the shock position is 
~ 0.025 eV/cm-' for 1 yuG magnetic field. This implies that, at 
the position of the shock, the main energy loss channel for rela- 
tivistic electrons is inverse Compton scattering in the CMB pho- 
ton field, characterised by an energy density of ~ 0.25 eV/cm^. 
Optical and infrared radiation from cluster galaxies is totally 
negligible. 

In Fig. Qwe plot the energy loss rates for a fiducial galaxy 
cluster. The thin lines represent the acceleration rates for a ve- 
locity Ml = 2000 km/s and different values of the magnetic 
field. The thick lines are the total energy loss rates for inverse 
Compton scattering off the CMB plus synchrotron emission. 
Calculations are done for three values of the upstream magnetic 
field Bi (1 yuG, 0.3 fiG, 0.1 yuG) and B2 = 4Bi. 

Below we consider the case of parameter values Bi - 0.3 /jG 
and Ml - 2000 km/s, which are reasonable for an accretion 
shock around a rich cluster (Blasi et al. 2007). In this case, the 
cut-off energy falls into the KN-dominated part of the total loss 
rate, at about 60 TeV (see Fig.|2ll. At the cut-off energy the loss 
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time scale is of the order of 4 x 10"^ yr. The evolution times 
of large-scale structures are several Gyr, meaning that electrons 
cool comparatively rapidly. Therefore the steady state is always 
reached. 

In Fig. [8^ we plot the electron spectrum at the shock, com- 
pared to the one expected for pure synchrotron losses. The mod- 
ification due to the dominant IC losses has the same characteris- 
tics as noticed before: a wider and less sharp cut-off. If we look 
at the effect this produces on the total spatially integrated parti- 
cle spectrum, we find that a small pile-up is still present around 
the cut-off energy (Fig. [SJ)). This feature is less pronounced in 
the present case because the energy losses right below the cut- 
off energy fall in the Thomson regime, as we can see from Fig.|2l 
This process does not happen deep enough in the Klein-Nishina 
regime to produce a strong feature. 

The photon spectra are affected in a similar way. In Figure 
we plot the upscattered spectrum of the CMB. We plot the 
two contributions to the total spectrum from the upstream and 
downstream regions, the latter dominating the total spectrum. 
Fig. |9}3 shows the synchrotron spectrum where the small pile-up 
feature of the electron spectrum is reproduced. 

4.3. Microquasars 

Another class of objects of interest are the microquasars. These 
are galactic binary systems composed of a regular O/B star being 
accreted onto a compact object (neutron star or black hole) that 
presents a jet. The nature of these objects is very complex and 
is not yet fully understood. Nevertheless the basic ingredients 
present in the system are a very intense stellar radiation field; 
e.g., for LS 5039, the luminosity of the star is L» ^ 7 x lO-'^ erg/s 
implying a radiation energy density up to Urad - 1000 erg/cm"*. 
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Fig. 11. a) Electron spectrum at the shock for a shock velocity 
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IC losses for an acceleration time > 100s, the dot-dashed one is 
for tacc - 50s, while the dashed one is obtained with pure syn- 
chrotron cooling in steady state, b) Integrated electron spectra 
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depending on the loc ation, at a temperature ofr^3.8xlO'*K 
(ICasares et al.ll2005l) . and a jet that is thought to be the site of 
particle acceleration. If blobs of plasma are emitted in the jet 
at slightly different velocity, internal sub-relativistic shocks can 
form, where electron acceleration can take place. The value of 
the magnetic field is not known, but has a reasonable range be- 
tween 0.01 and 1 G. 

For our calculation we refer to the model used in 
Khangul van et al.l (l2008l) , where the jet from the compact object 
is assumed to have an axis perpendicular to the orbital plane, 
and we put the acceleration zone at a distance Zq - 2x 10^^ cm 
from the jet base (corresponding to a distance equal to the mean 
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Fig. 12. Photon spectra produced by the electron distributions in 
Fig. [TTb . a) IC spectra for a diluted black body of temperature 
r = 3.8 X IC* K. b) Synchrotron spectra. 



orbital radius). At such a distance, the radiation energy density 
results: 



230 



erg 



(12) 



In Fig. [To] we plot the total energy loss rates due to IC plus 
synchrotron losses and the acceleration rates for a shock velocity 
of 30000 km/s, for three values of the upstream magnetic field 
Bi = 1 G, 0.3 G, 0.1 G, and a downstream magnetic field B2 - 
4Bi. We chose mi = 30000 km/s, corresponding to 10% of the 
speed of light, in order to remain consistent with the assumption 
of nonrelativistic shock. 

In this case the radiation energy density is considerably 
higher than the magnetic one, even for a field Bi = 1 G. This is 
the value we considered in our calculation in order to obtain an 
efficient acceleration of electrons up to TeV energies. Moreover, 
from a theoretical point of view, this choice allows us to explore 
a different region of the KN losses, compared with the previous 
two cases. 



The IC cooling time at the cut-off energy for the configura- 
tion described is - 40 s. As a rough estimate, we can say 
that the dynamical time of the shock is Tage ~ Zq/ui ^ 300 s. 
We therefore report on the results of our calculation for a time 
of 50 s, 100 s, and 300 s. We also show the steady state solu- 
tion for comparison. In Fig.[TT|we plot the electron spectrum at 
the shock and integrated over space. The difference in the shape 
of the spectrum around the cut-off energy between the expected 
spectrum for pure synchrotron cooling and the one obtained for 
dominant IC losses in panel a is very pronounced. We see that, 
when the age of the system is comparable to the loss time scale 
at the cut-off energy (50 s), the steady state is not reached even 
at the high energy end of the spectrum (Fig.[TT] both panels). In 
the cases corresponding to 100 s and 300 s, on the other hand, 
the spectrum around the cut-off reaches its steady state. Only the 
low energy part of the integrated spectrum presents the charac- 
teristic uncooled slope (Fig. [TTb). 

In this case the predominance of the IC is such that the hard- 
ening of the integrated spectrum is very pronounced, the pile-up 
spanning two orders of magnitude. 

As in the previous cases, the IC photon spectrum does not 
show prominent features (Fig. [T2h). while the synchrotron one 
reproduces the electron pile-up (Fig.[T2b). 



5. Conclusions 

We have addressed the problem of particle acceleration at non- 
relativistic shocks in the presence of strong energy losses, specif- 
ically in the never before treated case of electrons undergoing 
severe inverse Compton losses in the Klein-Nishina regime. 

We numerically solved the complete transport equation for 
the accelerated particles and calculated the particle distribu- 
tion function f{x, p) and the resulting radiation spectra. Our ap- 
proach provides a general tool for what concerns the species of 
particles accelerated (electrons and protons), the type of losses 
and the form of the diffusion coefficient. We focused on electrons 
undergoing synchrotron and inverse Compton losses, where the 
radiation energy density is higher than the magnetic one. We as- 
sumed Bohm diffusion: D(p) - pc^/(3eB). 

For high-energy electrons, with energy E > (mc^)^/eph, 
the inverse Compton scattering enters the Klein-Nishina regime 
where the cross-section decreases with energy. We evaluated the 
particle spectrum at the shock and found that, compared to the 
case of Thomson losses, the KN effect results in a broadening 
of the cut-off region and a hardening of the spectrum around 
the cut-off energy. This effect is most evident in the spatially in- 
tegrated spectrum, which exhibits a pronounced pile-up feature 
just below the cut-off energy. Such a feature can be significant, 
with an enhancement of the spectrum up to a few orders of mag- 
nitude, depending on the ratio of energy densities UradlU„,ag, the 
magnetic field strength, and the shock velocity. 

Once the electron spectrum is obtained, we can calculate the 
resulting spectra of the radiation, both synchrotron and inverse 
Compton. The hardening in the electron spectrum is barely vis- 
ible in the IC spectrum since the effect of the KN cross-section 
on the scattering process itself compensates the opposite action 
on the electron spectrum. The pile-up feature is reproduced by 
the synchrotron distribution, which no longer follows a simple 
power law below the cut-off, as in the case of dominant Thomson 
losses. 

We demonstrated the importance of this effect in three ex- 
amples of astrophysical objects in which the energy density of 
the radiation field dominates the magnetic field energy density. 
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namely a supernova remnant located in the Galactic centre re- 
gion, the accretion shocks of clusters of galaxies and the internal 
shocks in the jet of a microquasar, and discussed the possible 
consequences of our results on the emission of these objects. 
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